Phonon Softening near structural transition in BaFe 2 As 2 observed by inelastic x-ray 

scattering 
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In this work we present the results of an inelastic x-ray scattering experiment detailing the be- 
havior of the transverse acoustic [110] phonon in BaFe2As2 as a function of temperature. When 
cooling through the structural transition temperature, the transverse acoustic phonon energy is 
reduced from the value at room temperature, reaching a maximum shift near inelastic momentum 
transfer q = 0.1. This softening of the lattice results in a change of the symmetry from tetragonal 
to orthorhombic at the same temperature as the transition to long-range antiferromagnetic order. 
While the lattice distortion is minor, the anisotropy in the magnetic exchange constants in pnictide 
parent compounds is large. We suggest mechanisms of electron-phonon coupling to describe the 
interaction between the lattice softening and onset of magnetic ordering. 
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PACS numbers: 63.20.-e, 74.25.kc, 63.20.kd 



I. INTRODUCTION 

Discovery of superconductivity in iron pnictide 
materials 1 spurred intense experimental and theoretical 
effort to unravel the underlying physics of these com- 
pounds. The presence of FePn^ tetrahedra is universal 
in pnictide compounds, where iron atoms are arranged 
in a square planar lattice and tetrahedrally coordinated 
by pnictogen atoms. A common dynamic feature in pnic- 
tide families is a structural phase transition concomitant 
or in close vicinity with a transition to long range anti- 
ferromagnetic (AFM) order upon cooling 2 ^. This tran- 
sition is suppressed with hole^— ,electron2^, or isoelec- 
tronic doping^iS, eventually giving way to superconduc- 
tivity. 

At room temperature, the pnictide parent compound 
BaFe2As2 is a tetragonal (IA/mmm) paramagnet. At T 
sa 140 K 3 , BaFe2As2 undergoes a structural transition to 
orthorhombic (Fmmm) symmetry, and magnetically or- 
ders with spins AFM along the a-axis in the iron plane, 
ferromagnetic along the 6-axis, and AFM between the 
iron layers with an ordered moment of 0.87(3) (is per 
Fe siteii. The measured moment is lower than that pre- 
dicted by DFT calculations^. While the structural dis- 
tortion in the pnictides is small, the magnetic exchange 
anisotropy is very large when the exchange constants are 
calculated from an anisotropic Heisenberg model calcu- 
lation which includes damping with very small interlayer 
magnetic couplin g 13 ' 14 . 

For AFM ordering to occur, the high temperature 
tetragonal symmetry of BaFe2As2 must be broken to al- 
low for the spin anisotropy, making the lattice distortion 
from tetragonal to orthorhombic symmetry a require- 
ment for AFM ordering. The precise mechanism of the 
lattice distortion is a matter of considerable debate^— . 

Early theoretical work shows that the lattice 



distortion 2 in the oxyarsenide LaFcAsO is the manifesta- 
tion of relieved magnetic frustrations as the AFM order- 
ing results in different occupancy for the d xz , d yz orbitals 
which breaks the tetragonal symmetry^. However, the 
exchange constants determined from the spin wave dis- 
persion are not consistent with this idea^. Another ex- 
planation for the coupling of the transitional behaviors 
is that of ferro-orbital ordering. A minimal approach 
applying the Hartree-Fock approximation to a two or- 
bital model shows that AFM ordering arises as a conse- 
quence of ferro-orbital ordering, which then results in a 
lattice distortion^ 3 -. Three- and five-orbital models sug- 
gest that orbital degrees of freedom are strongly coupled 
to magnetic ordering, with d xz and d yz susceptible to 
orbital ordering, but that the orbital magnetization is 
much greater for the d yz orbital, resulting in a reduc- 
tion of the orbital hybridization between the d xz and d yz 
orbitals 2 ^. Other reports suggest that orbital ordering 
is able to spontaneously break the C4 symmetry locally 
without Fermi surface nesting or magnetic frustrations 21 . 

Softening of the shear elastic constant in 
Ba(Fe a: Coi_ a ;)2As2 has been noted in several works 
including resonant ultrasound spectroscopy on BaFe2As2 
and Ba(Feo.92Coo.o8)2As2r^ an d ultrasonic pulse echo 
measurements on undoped, underdoped 2 ^, and optimally 
doped Ba(Fe 2; Coi_ a ;)2As2 1 ^. The measurements showed 
a softening of the C^s mode near the structural or 
superconducting transition in all compounds where 
measurement wasn't precluded by sample instability 
against the mounting apparatus. The C*66 mode is 
responsible for structural phase transitions of type which 
serve to lower crystalline symmetry from tetragonal 
to orthorhombi o 27 ' 28 , so this observation is expected, 
although in the presented works several additional 
mechanisms for the phase transition are proffered. 

Resonant ultrasound spectroscopy measurements per- 
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formed on BaFe2As2 show a strong softening of the shear 
clastic constant, which is used to conclude that there is 
general lattice softening through the structural transi- 
tion driven by renormalization of the shear elastic con- 
stant by nematic fluctuations^ 5 -. In this context, the 
softening of the lattice is a secondary effect, with elec- 
tronic degrees of freedom taking precedence over elastic, 
and the structural transition then a mere consequence 
of the magnetic ordering. In the work conducted in ul- 
trasonic pulse echo spectroscopy, the softening minima 
tracks with the superconducting temperature under ap- 
plication of a magnetic fiel d ' 26 , and other elastic con- 
stants show almost imperceptible deviations from mono- 
tonic behavior in close proximity to the superconducting 
transition temperature 1 ^. One hypothesis derived from 
these experiments is the phase transition is driven by 
ferro-orbital ordering, suggestive of strong lattice cou- 
pling and orbital fluctuations 26 , while another concludes 
that the softening is induced by the coupling of the de- 
generate electronic orbitals to the elastic strain^. 

There is some disagreement between the reports on the 
clastic constant behavior at low temperature. Upon cool- 
ing below the structural or superconducting temperature, 
some reports show the value of Cqq recovering some 10% 
of the room temperature value at low temperature for the 
undoped and superconducting samples 25 , while other re- 
ports on the underdoped show no further observed soft- 
ening below the structural transition temperature^. 

Experimental studies on the transverse acoustic (TA) 
modes of the related compound CaFe2As2 shows a soften- 
ing of the transverse acoustic mode at the zone boundary 
for a TA [110] mode, and a mid-zone point at TA [100] 
at the magneto-structural transition when measuring ac- 
companied by phonon broadening 2 -. Further, calculation 
shows that phonons are very sensitive to the structural 
details of the FePn^ layers, with phonon modes that dis- 
tort the tetrahedra having substantial impact on the elec- 
tronic density of states close to Ep^. Anomalous phonon 
behavior is observed in Raman spectroscopy, with the E g 
in-plane phonon mode split in BaFc2As2^ r — , which is 
linked to strong spin-phonon coupling. This E g mode is 
a doubly-degenerate mode impacting Fe and As, and the 
contention is that the split is too large to be accounted 
for by the orthorhombic distortion and instead is indica- 
tive of spin-phonon coupling. Measurements of phonon 
behavior conducted with inelastic x-ray spectroscopy are 
replicated well by theory when magnetism is explicitly 
taken into account 54 , indicating the presence of spin- 
phonon coupling. 

The degree of hybridization between the As-4p and Fe- 
3d states is shown to have strong impacts on the Fe 2+ 
spin stated 5 -, and to influence the Fe-As phonon modest. 
Further, the dependence of the superconducting critical 
temperatur o 36 ' 37 and magnetic momen t 38 ' 39 of pnictide 
materials as a function of Fe-As layer separation is well 
established. 

The preceding discussion sets forth a clear need to un- 
derstand the specifics of the structural transition as they 



relate to the elastic degrees, magnetic, and orbital de- 
grees of freedom. Here we report the results of inelastic 
x-ray spectroscopy (IXS) measurements on single crys- 
tals of BaFe2As2 as a function of temperature. Since 
the acoustic phonon frequency goes as the square root 
of the elastic constant, we expect a softening of the low- 
energy acoustic phonon modes of BaFe2As2 when passing 
through the structural transition temperature. Phonon 
calculations predict a linear dispersion in the low-energy, 
long wavelength regim o 30 ' 34 . The acoustic phonon modes 
in BaFe2 As2 can be characterized with energy transfers of 
less than 10 meV, and the small and constant momentum 
resolution offered by IXS allows precise measurements of 
the phonon dispersion as a function of temperature. The 
previous measurements on the shear elastic constant 25 
provide information only at the zone center (r point, q 
= 0), so the IXS measurements described in this work 
provide additional information on the low-energy lattice 
dynamics in the BaFe2As2 system. 



II. EXPERIMENT 

We investigated the behavior of the transverse acous- 
tic phonon mode in BaFe2As2 near the (2,2,0)o position, 
which is the [110] mode in the reduced scheme. The 
orthorhombic notation and the Fmmm symmetry are 
used throughout this report. The measurements were 
made using the HERIX instrument on beamline 3-ID 
of the Advanced Photon Source at Argonne National 
Laboratory^ - 41 .. The HERIX instrument is on an inser- 
tion device beamline with a six-bounce monochromator 
setup providing an incident energy of 21.657 keV, and a 
resolution of 2.3 meV at the elastic line. 

Our sample was a single crystal of BaFe2 As2 measuring 
1mm x 2mm x 60 /im, grown using a self- flux method 7 . 
Calculated optimal thickness for a transmission mode ge- 
ometry for an incident energy of 21.657 keV was 55 /im. 
The crystal was oriented in transmission mode with (2, 
2, 0) and (1, 1, 3) in the scattering plane. This alignment 
allows measurement of the phonons in the basal plane of 
the form Q = G + q= (H-q,K + q), where G = (H, K) 
is the reciprocal lattice vector of the Bragg peak in the 
a — b plane, and q the inelastic momentum transfer. For 
temperature dependent studies, the crystal samples were 
mounted on copper holders with a small amount of var- 
nish, and the copper holders were attached to the cold 
finger of a closed cycle cryostat sealed by a beryllium 
dome. 

The beam size on sample was 50 ^m by 350 /im. Only 
one analyzer-detector pair was used for these measure- 
ments, and an analyzer mask was used to limit Bragg 
contamination for low q measurements. The mask size 
was 4 cm by 4 cm at q> 0.07, and 1 cm by 1 cm for all 
q < 0.07. The Stokes and anti-Stokes profiles were mea- 
sured in most cases, but at points where an anti-Stokes 
contribution to the phonon spectra was dynamically lim- 
ited, a contribution from elastic scattering was ensured 
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FIG. 1. (Color online) a) Comparison of 300 K, 140 K, and 
130 K phonon spectra at q — 0.1. b) Comparison of 300 K and 
130 K phonon spectra at q = 0.05. The dotted line shows the 
Bragg peak measurement at T=130 K, Q=(2,2,0)o- All spec- 
tra are Bose factor corrected and normalized to one. Errors 
are statistical errors propagated through Bose correction. 



by transversely scanning the crystal space at zero energy 
transfer. All data were normalized to the incident moni- 
tor count, a low efficiency ion chamber. 
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FIG. 2. (Color online) Energy of transverse acoustic phonon 
as a function of temperature for Q — (H — q,K + q) for 
H = K = 2 and q = (a) 0.15, (b) 0.12, (c) 0.10, (d) 0.08 in 
orthorhombic reciprocal lattice units. Dotted line indicates 
Ts- Error bars are statistical errors propagated through Bose 
correction and peak averaging. Two softenings are observed, 
one at Ts and additional softening at low q at low tempera- 
ture. 



For all temperatures, the phonon spectra were fit with 
Voigt profiles, employing the Whiting approximation to 
determine the widt h 42 i 43 . Phonon spectra were corrected 
for the Bose thermal population factors. 

Fitted phonon peak widths are constrained to that of 
a phonon measured at q = 0.2 at T = 130 K. Data for 
values of q > 0.07 were fit with two phonon peaks and 
an elastic contribution. The phonon energy for q > 0.07 
is determined by the average of the fitted peak positions 
of the Stokes and anti-Stokes contributions. 

Data for q < 0.07 at temperatures below 300 K can not 
be resolved as distinct peaks, so the width of the observed 
peak is used as a proxy for the degree of softening (Fig. 
[T]). The peak width is determined using a single Voigt 
profile. 



III. DISCUSSION 

The primary result of our investigation is the obser- 
vation of a reduction in energy of the [110]o transverse 
acoustic phonon in the low q region, evident from the 
measured phonon spectra (Fig. [T]). The data show this 
reduction is most pronounced at q = 0.1, reaching a max- 
imal amount near T = 140 K (Fig. As the tempera- 
ture is decreased, the phonon energy decreases, and then 
recovers, but not to the full energy observed at room 
temperature (Fig. [2} before decreasing again at low tem- 
peratures and low q. 

For data with q < 0.07 the fitted width of the peak 
centered around zero energy transfer exhibits a linear be- 
havior for 300 K (Fig. [3£i), with the width of the spectra 
collapsing into the Bragg peak width for lower tempera- 



4 




7^ 
6- 
5- 
4- 
3- 



I 2-.. 



0.00 



.J.-.!. 



• 


300K 


■ 


130K 


T 


57K 



0.02 0.04 0.06 0.08 0.10 
q(r.l.u-) 



FIG. 3. (Color online) a) Dispersion of TA phonon at 300 K, 
140 K, 130 K, and 57 K. Dotted line indicates a linear fit to 
the 300K data, b) Central peak width as a function of q for 
T=300 K, 130 K, and 57 K. Heavy dashed line at E = 2.3 
meV indicates the width of the Bragg peak measured at Q = 
(2,2,0) T=130 K. Dotted line on 300 K data represents the 
fitted peak width. Error bars are statistical errors propagated 
through Bose correction and the curve fitting procedure. 



tures (Fig. [3t>) . 

The lattice softening observed in the shear elastic mea- 
surements show a drop in the shear mode to a min- 
imal values at Tg, and then a slight recovery to low 
temperature o 16 ' 25 , or no recovery at all2S,. This is in con- 
trast to our observations here, as the phonon energy re- 
covers most of the energy after softening, and the slope 
of the dispersion changes at low temperature. 

The only mode investigated in this work was the low 
energy transverse acoustic mode, and there was only a 
local softening of the mode - the traditional linear dis- 
persion is recovered by q — 0.20 in all temperature cases 
(four temperature cases shown in Fig. [3]). This result 
is not captured in published DFT-GGA calculations re- 
viewed in preparation of this manuscript. 

The softening of the low energy transverse acoustic 
mode observed by IXS is expected in the context of soft 
mode transitions, as the decrease in energy of the long 
wavelength phonons is always seen in structural phase 
transitions of the type observed here 2 ^. But it has a few 
possible interpretations for microscopic mechanism. 

The first interpretation invokes Fermi surface nesting, 



which results in a regime of enhanced electron-phonon 
coupling. The Fermi surface structure of pnictides is well 
established to consist of hole pockets centered at the zone 
center, and electron pockets at the zone boundar y 44 ' 45 , 
with electron or hole doping serving to change the levels 
of the two pockets relative to one another—. The soften- 
ing of the TA phonon has a natural explanation in two- 
dimensional Fermi surface nesting giving rise to a Kolm 
anomaly, where strong quasiparticle excitations below 
q = 2kp result in enhanced electron-phonon coupling 47 . 
Indeed we observe softening from the zone center to a 
maximal phonon softening below q = 0.1, which is con- 
sistent with the value of 2kp = 0.1 for the inner hole 
pocket in BaFe2As2^- 

In a work on CaFe2As2, a series of line broadenings is 
observed 2 ^, though this is not used to draw the conclu- 
sion of Kohn anomalies, but may be related to the work 
here. However, possible broadening of the phonons for 
the BaFe2As2 compound studied here is too small to be 
fully resolved by the current technique. Line broadenings 
of this type have been observed in studies of elemental 
superconductors, where the phonon anomalies were fully 
resolved only when using high resolution spin-echo neu- 
tron spectroscopy 4 ^. In order to confirm the presence of 
phonon broadening and/or a Kohn anomaly in BaFe2As2, 
a careful analysis of the phonon lineshape using a higher 
energy resolution is required. Note also that there is some 
softening tendency at temperatures much below T s . The 
q range of this low temperature softening is distinct from 
the softening at T s ; below q = 0.12 at low temperatures, 
whereas softening is observed below q = 0.16 near T s . 
It appears that this softening has a different origin, and 
either one could be related to a Kohn anomaly. 

A second interpretation is to consider the phonon soft- 
ening in the context of coupling between the low energy 
acoustic mode and the orbital fluctuations. In this sce- 
nario, softening of the shear mode causes in-plane shifting 
of the iron atoms, impacting the interatomic distances 
along the ¥ePn^ layer. As these atoms shift their posi- 
tions, the tetragonal orbital degeneracy between the d xz 
and d yz orbitals is broken, allowing the manifestation of 
magnetic order. While the orthorhombic structural dis- 
tortion is small (of order 0.3%, with the change in the iron 
square lattice 2.8 A on a side to a rectangle 2.79/2.81 
A), and several works show that there is a strong de- 
pendence of the electronic structure on small details of 
the lattice, in particular the height of the pnictogen atom 
above the iron plane^ 5 - 2 . Additionally, it has been shown 
that there is an enhancement of electron-boson coupling 
arising from orbital anisotropy^i, behavior which is in 
contrast to lower mean field values for the estimate of 
such couplings 50,51 . Further, the pnictide systems are 
very close to the Stoner criticality, where even a small 
lattice distortion can have a large consequence to the 
magnetism in the system^. 

The proposed mechanism indicates that other phonon 
modes may be involved in the structural transition via 
a similar softening behavior^. Modes that directly in- 
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fluence the dynamics of the FeAs4 tetrahedra are being 
examined. Further, a coupling of acoustic and optical 
modes as is seen in other structural phase transitions 53 
is not ruled out, and requires further examination. 



fully. This observation is explained in terms of electron- 
phonon coupling resulting in a lifting of the orbital de- 
generacy in the paramagnetic tetragonal state, allowing 
full manifestation of the AFM order. 



IV. CONCLUSIONS 

We have investigated the behavior of the transverse 
acoustic phonon in BaFe2As2 using IXS, and found that 
this phonon mode softens as a function of temperature, 
with a maximal softening near q = 0.1. The energy of 
the phonon mode drops well below that observed at room 
temperature, reaching a minimum near the structural 
phase transition temperature before recovering almost 
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